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This work presents the study of the excited states of the Cooper problem in the three- 
dimensional Anderson model. It is shown that the excited pair states remain localized while 
their excitation energy AE is negative. For AE > the particles are delocalized over the 
' three-dimensional lattice. 
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The Cooper problemEl is the cornerstone of the well known BCS theory for superconductivity. 
Even if the Cooper problem deals with only two interacting particles (TIP) above a frozen 
Fermi sea, it nevertheless captures the essential physical features of the superconducting states. 
^ . Indeed, in comparison with the BCS theory, the solution of the Cooper problem leads to the 

appearance of coupled states with a qualitatively correct coupling energy and pair size. Then it 
seems natural that the study of the Cooper problem with disorder could provide a useful step 
• in the understanding of the superconductivity in presence of relatively strong disorder. The 

first studies s>n the Cooper problem with disorder were done for three-dimensional i and two 
' dimensional B systems. These studies era were focused on the properties of the ground state of 

■ two particles above a frozen Fermi sea coupled via an on-site attractive Hubbard interaction. 

This interaction creates a phase of bi-particle localized states (BLS) in the regime of disorder 
where non-interacting states are delocalized. At the same time the mean- field solution of the 
Cooper problem (Cooper ansatz) gives delocalized pairs. This shows that the non-diagonal 
matrix elements of interaction play an important role in presence of disorder. These TIP results 
are in qualitative agreement with recent many-body investigations §'! of the ground state of the 
attractive Hubbard model with disorder. 



Here the studies are concentrated on the properties of excited TIP states. For that purpose 



we introduce the Anderson Hamiltonian for one-particle, 



n) (m| (1) 

< I ■ n (n,m) 



where n and m are the index vectors on the three-dimensional lattice with periodic conditions, 
(, ) denotes nearest neighbor sites, V is the hopping term and the random on-site energies E n are 
homogeneoulsy distributed in the energy interval — ^- , where W is the disorder strength. 
To study the effects of the attractive interaction (U < 0) on two particles near the Fermi sea we 
generalize the Cooper approaclU for the disordered case. We write the TIP Hamiltonian in the 
basis of one-particle eigenstates of the Hamiltonian (|l|). In this basis the Schrodinger equation 
for TIP reads 

(E mi + E m2 ) X mi)1n2 + U ^ Qm 1 ,m 2 ,m' l ,m 2 X m \,m' 2 = ^\ Xmi.mj' (^) 

Here E m are the one-particle eigenenergies corresponding to the one-particle eigenstates \<j> m ) 
of Hi and Xmi m 2 are ^ ne components in the non- interacting eigenbasis | <p mi <g) 4> m2 ) of the A th 
TIP eigenstate x X ^ corresponding to the eigenenergy E\. The matrix elements 1) r Qm 1 ,m 2 ,m,' 1 .m, , 2 
give the interaction induced transitions between non-interactive eigenstates \(f) mi <g> (j>m 2 ) an d 
| 0m'j ® ^ • These matrix elements are obtained by writing the Hubbard interaction in the 
non-interactive eigenbasis of the model (|l]). In analogy with the original Cooper problemi the 



summation in @ is done over the non-interacting states above the Fermi level, in the labelling 
E m > > Ep corresponds to m'x 2 > 0. The Fermi energy is determined by a fixed filling factor 
v = 1/2. To keep similarity with the Cooper problem we restrict the summation on m' l 2 by the 
condition 1 < m! x + m' 2 < M. In this way the cut-off of M unperturbed orbitals introduces an 
effective phonon frequency u>o oc M/L 3 = 1/a where L is the linear system size. When varying 
L we keep a fixed so that the phonon frequency is independent of the system size. All the data 
in this work are obtained with a = 30 but we also checked that the results are not sensitive to 
the change of a. 

In order to analyze the properties of the TIP excited states of the disordered Cooper problem 
we compute the energy level spacing distribution P(s) obtained by the diagonalization of 
From -P(s) we calculate 

f°° [P( S )-P WD ( S )}ds 

Here Pp(s) = exp(— s) and Pwd(s) = ^exp^— ^j-^j are respectively the Poisson and the 
Wigner-Dyson distributions, and so — 0.4729 is their intersection point. In this way, r/ varies 
from 1 [P(s) = Pp(s)] to [P(s) = Pwd(s)] and thus characterizes a transition from localized 
to delocalized states. For example^ for the one-particle problem this method allows to detect 
efficiently the Anderson transition 
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Figure 1: Dependence of r\ on the excitation energy AE for (a) W = 18V > W c — 16.5V and 
(b) W = 12V < W c . The interaction strength is U = — 4V and the linear lattice size is L — 
6 (o), 7 (•), 8 (A), 10 (full triangle), 12 (square) and 14 (full square). The number of disorder realization used is 

N D = 3000 and a = 30. 

In Fig. g, r] is shown for disorder strength W = 18V (a) and W = 12V (b), an interaction 
strength U = —AV and for different linear lattice size, from L = 6 up to L = 14. 77 is presented 
versus the excitation energy AE = E — 2Ep where E is the energy of the two interacting 
particles. Without interaction for W = 18V > W c ~ 16.51/ ]] the one particle states are well 
localized (see Fig. |2|a). When the on-site Hubbard attraction is switched on, coupled states 
appear below the Fermi energy of two non-interacting particles (AE < 0). Fig. |l|a clearly shows 
that at the thermodynamic limit the coupled states with AE < are localized (rj — ► 1). As 



a At the band center W c — 16.5V is the critical value of the disorder at which the Anderson transition occur; 




W > W c all the non- interacting orbitals are localized in the lattice basis. The on-site attractive 
interaction acts as an additional constraint forcing the two particles (with opposite spin) to stay 
together in the same well of potential (see Fig. ^o). Above the Fermi energy, 77 tends slowly to 
indicating that these unbounded states (AE > 0) are delocalized. This derealization is due to 
the fact that an interaction (repulsive or attractive) between particles destroies single particle 
localization and leads to a propagation of pairs of size l\ over a distance l c much larger than l\ 
Utl. The enhancement factor k is then determined by the density P2 of two-particle states coupled 
by the interaction and the interaction induced transition rate T2 between noninteracting states, 
so that k = T2P2- For AE > the density p2 grows with the excitation energy (p2 °< AEp that 
strongly enhances the derealization of pairs (see Fig. ^c). 

In Fig. |]b, r\ is shown for W = 12V < W c . For this regime the non-interacting one-particle 
states are delocalized (see Fig. |2|d). For the different linear lattice size L the ground state 
is still localized (77 ~ 1) in agreement with!. Indeed as W/W c ~ 0.72 > W s /W c ~ 0.35 the 
TIP ground state is in the BLS phase, where W s is the critical value of the disorder for the 
superconductor-insulator transition found in@ for TIP pairs. In the region — 0.25V < AE < 
although 77 drops from 0.9 to 0.25 a certain localization of TIP pairs still remains (see Fig. |2|s). 
For AE > 0, the TIP states become delocalized (see Fig. ||f). 
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Figure 2: One-particle probability distribution /a(iii) = ^ n | Xni,n 2 | extracted from TIP eigenstates | \ X } 
and projected on a (x,y)-plane. The linear lattice size is L — 16. f\(n\) is presented for W = 18V [ (a) ground 
state without interaction, (b) excited states A = 25 with AE\ ~ —0.25V for U = ~4V and (c) excited states 
A = 758 with AE\ ~ 0.25V for U = -W } and for W = 12V [ (d) ground state without interaction, (e) excited 
states A = 83 with AE X ~ -0.05V for U = -IV and (f) excited states A = 512 with AE X ~ 0.125V for U = -4V]. 



In order to characterize the wave function properties of the excited states of the disordered 
Cooper problem, we compute 
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(4) 



E x &[E-e,E+e] 



where the brackets mark the averaging over = 400 realizations of the disorder and over an 
energy interval 2e centered on E. The inverse participation ratio (IPR) £ (E) counts the average 
number of non- interacting states \<f> mi <8> 4>m 2 ) occupied by the TIP wave functions belonging to 
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Figure 3: IPR £ as a function of the coupling energy AE for a disorder strength W — 12 (o) [A.E/V is then read 
on the upper horizontal axis] and W = 18 (•) [A.E/V is then read on the lower horizontal axis]. The interaction 
strength is U — —4V and the linear lattice size L — 12. Nd ~ 300 disorder realization are used, a = 30, 2e ~ 0.05 

for W = 18 and 2e ~ 0.025 for W = 12V. 

energy interval [E — e, E + e]. Fig. || shows the IPR as a function of the coupling energy AE 
for W > W c and W < W c . For W = 18V as the one-particle noninteracting states are localized, 
only few of these states (£ ~ 5 — 10) are enough to build the TIP localized states with AE < 0. 
On the contrary for W = 12V the TIP localized states (AE < 0) occupy more noninteracting 
states (£ ~ 30) than for W = 18V (Fig. ||). Indeed, for W < W c the one-particle noninteracting 
states are delocalized and it is necessary to have a sufficient number of these states to see the 
TIP eigenstates localization induced by the quantum interferences. 

In conclusion, the studx_pf the excited states of the disordered Cooper problem shows that 
the BLS phase found inH'Ba exists for TIP eigenstates with energy AE < 0. This phase is 
destroyed for TIP eigenstates above the Fermi level AE > 0. 
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